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Abstract Naively, the "best" method of renormalization is the one where a momentum cutoff is taken 
to infinity while maintaining stable results due to a cutoff-dependent adjustment of counterterms. 
We have applied this renormalization method in the non-perturbative calculation of phase-shifts for 
nucleon-nucleon (NN) scattering using chiral NN potentials up to next-to-next-to-next-to-leading 
order (N 3 LO). For lower partial waves, we find that there is either no convergence with increasing 
order or, if convergence occurs, the results do not always converge to the empirical values. For higher 
partial waves, we always observe convergence to the empirical phase shifts (except for the 3 Gs state). 
Furthermore, no matter what the order is, one can use only one or no counterterm per partial wave, 
creating a rather erratic scheme of power counting that does not allow for a systematic order-by-order 
improvement of the predictions. The conclusion is that infinite-cutoff renormalization is inappropriate 
for chiral NN interactions, which should not come as a surprise, since the chiral effective field theory, 
these interactions are based upon, is designed for momenta below the chiral-symmetry breaking scale 
of about 1 GeV. Therefore, this value for the hard scale should also be perceived as the appropriate 
upper limit for the momentum cutoff. 

Keywords Chiral perturbation theory • Nucleon-nucleon scattering • Non-perturbative renormaliza- 
tion 



1 Introduction 

During the past two decades, it has been demonstrated that chiral effective field theory (chiral EFT) 
represents a powerful tool to deal with hadronic interactions at low energy in a systematic and model- 
independent way (see Refs. [T|[2] for recent reviews). For the construction of an EFT, it is crucial to 
identify a separation of scales. In the hadron spectrum, a large gap between the masses of the pions 
and the masses of the vector mesons, like p(770) and a; (782), can clearly be identified. Thus, it is 
natural to assume that the pion mass sets the soft scale, Q ~ m^, and the rho mass the hard scale, 
A x ~ m p ~ 1 GeV, also known as the chiral-symmetry breaking scale. This is suggestive of considering 
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a low-energy expansion arranged in terms of the soft scale over the hard scale, (Q/ A X Y , where Q is 
generic for an external momentum (nucleon three-momentum or pion four-momentum) or a pion mass. 
The appropriate degrees of freedom are, obviously, pions and nucleons, and not quarks and gluons. 
To make sure that this EFT is not just another phenomenology, it must have a firm link with QCD. 
The link is established by having the EFT observe all relevant symmetries of the underlying theory, in 
particular, the broken chiral symmetry of low-energy QCD [3]. 

The early applications of chiral perturbation theory (ChPT) focused on systems like irir [3] and 
ttN [5], where the Goldstone-boson character of the pion guarantees that a perturbative expansion 
exists. But the past 15 years have also seen great progress in applying ChPT to nuclear forces [TJ [5J 
El [7J El [91 [101 HIl H2I HH [El HSl H2 [HI HH EUl EH EH E3] . The nucleon-nucleon (NN) system is 
characterized by large scattering lengths and bound states indicating the non-perturbative character 
of the problem. Weinberg [5J [7] therefore suggested to calculate the nuclear amplitude in two steps. In 
step one, the nuclear potential, V, is defined as the sum of irreducible diagrams, which are evaluated 
perturbatively up to the given order. Then in step two, this potential is iterated to all order (i.e., 
summed up non-perturbatively) in the Schrodinger or Lippmann-Schwinger (LS) equation: 

r(p p) = v( P ', P ) + / dV v(p', P ") M ^ 2 , f ( P ", P ) , (i) 

J p A — p" + ie 

where T denotes the NN T-matrix and Mjv the nucleon mass. 

In general, the integral in the LS equation is divergent and needs to be regularized. One way to do 
this is by multiplying V with a regulator function 

V( P p) —> V(p ', p) e-(*"M) 2 " e -(pM) 2 " . (2 ) 

Typical choices for the cutoff parameter A that appears in the regulator are A ~ 0.5 GeV <C A x ss 1 
GeV. 

It is pretty obvious that results for the T-matrix may depend sensitively on the regulator and its 
cutoff parameter. This is acceptable if one wishes to build models. For example, the meson models of 
the past Ell [55] always depended sensitively on the choices for the cutoff parameters which, in fact, 
were important for the fit of the NN data. However, the EFT approach wishes to be fundamental in 
nature and not just another model. 

In field theories, divergent integrals are not uncommon and methods have been developed for how 
to deal with them. One regulates the integrals and then removes the dependence on the regularization 
parameters (scales, cutoffs) by renormalization. In the end, the theory and its predictions do not depend 
on cutoffs or renormalization scales. So-called renormalizable quantum field theories, like QED, have 
essentially one set of prescriptions that takes care of renormalization through all orders. In contrast, 
EFTs are renormalized order by order. 

The renormalization of perturbative EFT calculations is not a problem. The problem is non- 
perturbative renormalization. This problem typically occurs in nuclear EFT because nuclear physics is 
characterized by bound states which are non-perturbative in nature. EFT power counting may be dif- 
ferent for non-perturbative processes as compared to perturbative ones. Such difference may be caused 
by the infrared enhancement of the reducible diagrams generated in the LS equation. 

Weinberg's implicit assumption [51 [TJ [55] was that the counterterms introduced to renormalize the 
perturbatively calculated potential, based upon naive dimensional analysis ("Weinberg counting"), are 
also sufficient to renormalize the non-perturbative resummation of the potential in the LS equation. In 
1996, Kaplan, Savage, and Wise (KSW) [5711551153] pointed out that there are problems with the Wein- 
berg scheme if the LS equation is renormalized by minimally-subtracted dimensional regularization. 
This criticism resulted in a flurry of publications on the renormalization of the non-perturbative NN 
problem [I2[31][32[33l[10[33[^ 

The literature is too comprehensive to discuss all contributions. Let us just mention some of the work 
that has particular relevance for our present discussion. 

If the potential V consists of contact terms only (a.k.a. pion-less theory), then the non-perturbative 
summation Eq. ([I]) can be performed analytically and the power counting is explicit. However, when 
pion exchange is included, then Eq. ([I]) can be solved only numerically and the power counting is 
less transparent. Perturbative ladder diagrams of arbitrarily high order, where the rungs of the ladder 
represent a potential made up from irreducible pion exchange, suggest that an infinite number of 
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counterterms is needed to achieve cutoff independence for all the terms of increasing order generated by 
the iterations. For that reason, KSW [27| [28J [29] proposed to sum the leading-order contact interaction 
to all orders (analytically) and to add higher-order contacts and pion exchange perturbatively up to 
the given order. Unfortunately, it turned out that the order by order convergence of one-pion exchange 
(1PE) is poor in the 3 Si- 3 Dx state [5D1ISI]- The failure was triggered by the 1/r 3 singularity of the 
1PE tensor force when iterated to second order. Therefore, KSW counting is no longer taken into 
consideration (see, however, [27]). A balanced discussion of possible solutions can be found in [53] . 

Some researchers decided to take a second look at Weinberg's original proposal. A systematic 
investigation of Weinberg counting in leading order (LO) has been conducted by Nogga, Timmermans, 
and van Kolck [3DJ in momentum space, and by Valderrama and Arriola at LO and higher orders in 
configuration space [351 1371 39 . A comprehensive discussion of both approaches and their equivalence 
can be found in [131 HE] • 

The LO NN potential consists of 1PE plus two nonderivative contact terms that contribute only 
in S waves. Nogga et al find that the given counterterms renormalize the S waves (i.e., stable results 
are obtained for A — > oo) and the naively expected infinite number of counterterms is not needed. 
This means that Weinberg power counting does actually work in S waves at LO (ignoring the m r 
dependence of the contact interaction discussed in Refs. [27] [25] [2D] E3] ) • However, there are problems 
with a particular class of higher partial waves, namely those in which the tensor force from 1PE is 
attractive. The first few cases of this kind of low angular momentum are 3 Pq, 3 P2, and 3 I?2, which 
need a counterterm for cutoff independence. The leading order (nonderivative) counterterms do not 
contribute in P and higher waves, which is why Weinberg counting fails in these cases. But the second 
order contact potential provides counterterms for P waves. Therefore, the promotion of, particularly, 
the 3 Pq and 3 P2 contacts from next-to-leading order (NLO) to LO would fix the problem in P waves. 
To take care of the 3 Z>2 problem, a next-to-next-to-next-to-leading order (N 3 LO) contact needs to be 
promoted to LO. Partial waves with orbital angular momentum L > 3 may be calculated in Born 
approximation with sufficient accuracy and, therefore, do not pose renormalization problems. In this 
way, one arrives at a scheme of 'modified Weinberg counting' |36j for the leading order two-nucleon 
interaction. 

For a quantitative chiral NN potential one needs to advance all the way to N 3 LO [2D]. Thus, 
the renormalization issue needs to be discussed beyond LO. Naively, the most perfect renormalization 
procedure is the one where the cutoff parameter A is carried to infinity while stable results are main- 
tained. This was done successfully at LO in the work by Nogga et al [3DJ described above. At NNLO, 
the infinite-cutoff renormalization procedure has been investigated in [43; 44; 45j for partial waves with 
total angular momentum J < 1 and in [39] for all partial waves with J < 5. At N 3 LO, only a study of 
the state exists [32]. Thus, a full analysis of the issue is still lacking. 

It is, therefore, the purpose of this paper to study the method of (non-perturbative) infinite-cutoff 
renormalization systematically order-by-order from LO to N 3 LO and for all partial-waves with J < 6. 
As discussed, it is necessary to carry this investigation through all of these orders, because the presently 
existing quantitative chiral NN potentials are of order N 3 LO and their renormalizability needs to be 
investigated. 

This paper is organized as follows. In Sec. 2, we will present the chiral NN potential up to order 
N 3 LO and, in Sec. 3, the non-perturbative renormalization will be discussed. The order by order 
convergence (or non-convergence) is the subject of Sec. 4, and Sec. 5 will conclude the paper. 



2 The chiral NN potential up to N 3 LO 

EFTs are defined in terms of effective Langrangians which are given by an infinite series of terms 
with increasing number of derivatives and/or nucleon fields, with the dependence of each term on the 
pion field prescribed by the rules of broken chiral symmetry. Applying this Lagrangian to a particular 
process, an unlimited number of Feynman graphs can be generated. Therefore, we need a scheme that 
makes the theory manageable and calculabel. This scheme which tells us how to distinguish between 
large (important) and small (unimportant) contributions is ChPT, and determining the power v of the 
expansion has become known as power counting. 
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Nuclear potentials are denned as sets of irreducible graphs up to a given order. The power v of a 
few-nucleon diagram involving A nucleons is given in terms of naive dimensional analysis by: 

v = -2 + 2A-2C + 2L + (3) 



ith 



11 ■ 

A = d l + ^~2, (4) 



where C denotes the number of separately connected pieces and L the number of loops in the diagram; 
di is the number of derivatives or pion-mass insertions and rn the number of nucleon fields (nucleon 
legs) involved in vertex i; the sum runs over all vertices contained in the diagram under consideration. 
Note that Ai > for all interactions allowed by chiral symmetry. For an irreducible NN diagram 
( "two-nucleon potential", A = 2, C = 1), Eq. ^ collapses to 

v = 2L + J2 A <- ( 5 ) 

i 

Thus, in terms of naive dimensional analysis or "Weinberg counting" , the various orders of the 
irreducible graphs which define the chiral NN potential are given by: 

Vio = Kf + ^ 0) (6) 
V NLO = V w + V™ + V^> + V£> (7) 

V N NLO = V N LO + V$ + V$ (8) 
V N 3 LO = V NNLO + Kt 4 ' + V$ + + V£ ] (9) 

where the superscript denotes the order v of the low-momentum expansion. LO stands for leading order, 
NLO for next-to- leading order, etc.. Contact potentials carry the subscript "ct" and pion-exchange 
potentials can be identified by an obvious subscript. 
The charge-independent 1PE potential reads 

Vi*(p ,P) = —rk T i ' T 2 y~ — 2 — > ( 10 ) 

where p ' and p designate the final and initial nucleon momenta in the center-of-mass system (CMS) 
and q = p ' — p is the momentum transfer; (T1.2 and are the spin and isospin operators of nucleon 
1 and 2; gA, f v , and m T denote axial- vector coupling constant, the pion decay constant, and the pion 
mass, respectively. We use f„ = 92.4 MeV and gA = 1-29 throughout this work. Since higher order 
corrections contribute only to mass and coupling constant renormalizations and since, on shell, there 
are no relativistic corrections, the on-shell 1PE has the form Eq. ( JTo] ) in all orders. 

In this paper, we will specifically calculate neutron-proton (np) scattering and take the charge- 
dependence (isospin violation) of the 1PE into account. Thus, the 1PE potential that we actually 
apply reads 

V^ p) (p ',p) = -V l7r (rM + (-l) T+1 2V llT ( m7I ±) , (11) 
where T denotes the isospin of the two-nucleon system and 

Vi„{m n ) = - 772- — 2^ 2~ ■ ( 12 ) 

4/2 q2 + m 2 



We use m,o = 134.9766 MeV and m,± = 139.5702 MeV. 
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2.1 Leading order (LO) 

The LO chiral NN potential consists of a contact part and an 1PE part, cf. Eq. The 1PE part is 
given by Eq. (11) and the LO contacts are 

y c ( t 0) (p',p) = C s + C T( T 1 . ( T 2 , (13) 

and, in terms of partial waves, 

V^\ 1 S )^C 1So ^4n(C s -3C T ) 

K ( t 0) ( 3 ^) = Cs Sl = 4tt (C s + C T ) , (14) 
where Cs, Ct, Cis > C 3 Si are constants. 



2.2 Next-to-leading order (NLO) 

Multi-pion exchange starts at NLO and continues through all higher orders. It involves divergent loop 
integrals that need to be regularized. An elegant way to do this is dimensional regularization which 
(besides the main nonpolynomial result) typically generates polynomial terms with coefficients that 
are, in part, infinite or scale dependent [13J . One purpose of the contacts is to absorb all infinities 
and scale dependencies and make sure that the final result is finite and scale independent. This is the 
renormalization of the perturbatively calculated NN potential, which must be carefully distinguished 
from the non-perturbative renormalization to be discussed in Sec. [3j The perturbative renormalization 
of the NN potential is very similar to what is done in the ChPT calculations of irir and irN scattering, 
namely, a renormalization order by order, which is the method of choice for any EFT. 

For the NLO chiral NN potential, Eq. we need to specify the second order contact part and 
the two-pion exchange (2PE) part. The NLO contact terms are given by [T] 

K ( t 2 V,p) = c ig 2 + c 2 k 2 

+ (C 3 q 2 + C A k 2 )<r 1 -* 2 
+ C 5 (-iS • (q x k)) 
+ C 6 (<ti • q) (er 2 • q) 

+ C 7 K ■ k) (<r 2 • k) , (15) 

with the partial-wave decomposition 













= C 3 p pp' 






= Ci Pl pp' 






= Cs Pl pp' 






= C 3Sl (p 2 + P ' 2 




( 3 ^l- 3 ^l) 


= C3 Sl ^3 Dl p 2 




( 3 £>i- 3 Si) 


= C3 Sl -3 Dl P 






= C 3 p 2 pp' . 



To state the 2PE potentials, we introduce the following notation: 
V 2n (p ',p) - V c +t 1 -t 2 W c 

+ [Vs + Ti • T 2 Ws ] <Ti ■ <T 2 

+ [V LS + ti • t 2 W LS ] (-iS • (q x k) ) 
+ [ Vt + i~i • r 2 Wt ] ci • q cr 2 • q 

+ [V a L + Ti ■T 2 W aL } o-i-(qxk) (T 2 -(qxk), (17) 
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where 



k = i (p ' + p) is the average momentum, and 

S = ±(«7i 



2 V" i i C2) the total spin. 
Using the above notation, the NLO 2PE is simply given by [TJ [T3 



where 



and 



Wg>=- 



v, 



(2) 



£(g) 

3847r 2 /^ 



4m 2 (5 ff l - 4 5 2 - 1) + ^(23^ - 10<£ - 1) 
64tt2/4 ' 



— In 



q 2m,r 
= y/Aml + q 2 . 



(18) 

(19) 
(20) 

(21) 
(22) 



2.3 Next-to-next-to-leading order (NNLO) 



There are no new contacts at NNLO, cf. Eq. pi), and, thus, all we need is the third order 2PE potential, 
which is [using the notation introduced in Eq. (17l] [T1IT3] 



where 



and 



.(3) 
"C 



w, 
v. 



(3) 

ci 

(3) 
Tl 



,(3) 
"C2 ) 

w (3) 

.(3) 
"T2 ' 



Tl 



(3) 



T2 ' 



w. 



r 3) = 

j. 3) = M 

(3) 

(3) _ w (3) 

S ~ yy SI ^ " S2 ' 



y(3) , y (3) 
V S1 + V S2 ) 



(3) 



(3) = 3g> 2 A(g) 
LS 32irM N f* 

■•1 _ Sa(I-Sa) 



^ 2 A(g) , 



y(3) 



16tt/. 



32nM N f* 



2mi(2cx-c 3 )-q 2 c 3 + 



35i 



w 2 A( g ) 



1287^/4 



{3 5 ^m> 2 - [4m 2 



16Mjv 

2g 2 - <&(4m» + 3q 2 )} w 2 A(q)} , 



1 T/(3) _ 9g\w 2 A(g) 



-V y 



ir. 



(3) _ 



512irM N f± ' 
32^/4 



C4 



1 

4Mjv 



^(lOrn 2 



3<? 2 ) 



V C2 — 



35l 



(3) 



C2 



T/ (3) - 

Vrprf 



256tt/4M w 



1 



-v {3) = 
f Vs2 5127T/4M 



'T2 



IK 



(3) 



(3) 



T2 



-2^S2 



(m^w 2 +w 4 A(q)), 
(m^w 2 + w 4 A(q)) , 

(m„ + w 2 A(q)) , 
K + ™ 2 -4(<?)) 



35l 



N 



256tt^M n 



(23) 
(24) 
(25) 
(26) 
(27) 
(28) 

(29) 
(30) 

(31) 
(32) 
(33) 
(34) 

(35) 
(36) 
(37) 
(38) 
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with 



and 



A(q) = — arctan - - 
w 2q 2m v 



'2ml + 1 



(39) 
(40) 



Equations (35 1- (38 1 are corrections of the iterative 2PE, see Ref. [T] for details. In all 2PE potentials, we 
use the average nucleon mass, Mm = 938.9182 MeV, and the average pion mass, — 138.039 MeV. 
The values for the low-energy constants are c\ = —0.81 GeV -1 , C3 = —3.20 GeV -1 , and C4 = 5.40 
GeV" 1 pp. 



2.4 Next-to-next-to-next-to-leading order (N 3 LO) 
At N 3 LO, 14 new contact terms appear pQ, 

Kt 4) &>', P) = D iq A + D 2 k 4 + D 3 q 2 k 2 +U 4 (qx k) 2 

+ (D 5 q 4 + D 6 k 4 + D 7 q 2 k 2 + D 8 (q X k) 2 ) a x • er 2 
+ (D 9 q 2 + D 10 k 2 ) HS-(qxk)) 
+ (flu q 2 + D12 k 2 ) (<ri • q) (<r 2 ■ q) 
+ (D la q 2 + £ 14 fc 2 ) (c^ • k) (o-2-k) 

+ Di8 (o-i • (q X k) <r 2 • (q X k)) , (41) 
which contribute as follows to the partial- wave potentials, 



VP?*) 


= D lSa (p' 4 +p 4 ) + D 1 


& 2 

SoP P 




= D 3Pa (p' 3 p + p'p 3 ) 






= D 1Pl (p' 3 p + p'p 3 ) 






= D 3Pl (p' 3 p + p'p 3 ) 






= D 3Sl (p' 4 +p 4 ) + D 3 


/2 2 
SiP P 




= D 3Dl p' 2 p 2 






= D 3Sl _3 Dl p 4 +Ds Sl . 


/2 2 
-*DiP P 


vPflh-'Si) 


= Ds Sl _3 Dl p' + Ds Sl 


/2 : 

-*D t P P 


v c [ 4) Cd 2 ) 


= D lD2 p' 2 p 2 






= D 3D2 p' 2 p 2 






= D SP2 (p' 3 p + p'p 3 ) 




V^( 3 P 2 -*F 2 ) 


= D 3 p 2 _ 3F2 p'p 3 




v£\ 3 F 2 -*P 2 ) 


= D 3P2 _ 3F2 p' 3 p 






= D 3D y 2 P 2 . 





The 2PE contributions at this order, V 2 \ , are very involved, which is why we will not reprint them 
here. The comprehensive expressions can be found in Appendix D of Ref. pQ. We note that, in the 
calculations of this paper, we apply all N 3 LO 2PE terms including those which require numerical 
integrations. The parameters we use in this work are listed in column "NN potential" of Table 2 of 
Ref. p. 

The N 3 LO three-pion exchange (3PE) contributions v£*\ cf. Eq. (El), are left out, since they have 
been found to be negligible [S3J |S3] ■ 
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3 NN scattering and non-perturbative renormalization 

For the unitarizing scattering equation, we choose the relativistic three-dimensional equation proposed 
by Blankenbecler and Sugar (BbS) [57], which reads, 

f r1 3 r>" M 2 1 

T(p',p) = V(p',p) + / ±P V (p',p")^ - „ r(p",p) (43) 

with E p ii = J Mjy + p" 2 . The advantage of using a relativistic scattering equation is that it automat- 
ically includes relativistic corrections to all orders. Thus, in the scattering equation, no propagator 
modifications are necessary when raising the order to which the calculation is conducted. 
Defining 

nP'»P)-7i^np',p)^ (44) 



(2tt) 3 V E p , ^ ' tr ' V E p 



and 



M N 



n "'^ s ^?\jt T(p '- p) it' (45) 

where the factor l/(27r) 3 is added for convenience, the BbS equation collapses into the usual, nonrel- 
ativistic Lippmann-Schwinger (LS) equation, Eq. ([l]). Since V satisfies Eq. ([lj, it can be used like a 
usual nonrelativistic potential, and T may be perceived as the conventional nonrelativistic T-matrix. 



The square-root factors in Eqs. (|44j) and (45) are applied to the potentials of all orders except in LO. 
In the LS equation, Eq. 0, we use 

2M v M n , . 

M N = — p " = 938.9182 MeV, and (46 
M p + M n v ' 

2 _ M p 2 T lab (71 ab + 2M») 

P ~ (M p + M„) 2 + 2T lab M p ' ( ' 

where M p = 938.2720 MeV and M n = 939.5653 MeV are the proton and neutron masses, respectively, 
and Ti a b is the kinetic energy of the incident neutron in the laboratory system ("Lab. Energy"). The 
relationship between p 2 and Xi a b is based upon relativistic kinematics. 

We renormalize the LO chiral NN potential as described in Refs. [3S1 |3H] and discussed in the 
Introduction. We then perform the infinite-cutoff renormalization also for the NLO, NNLO, and N 3 LO 
chiral NN potentials. This is accomplished by studying the dependence of the phase shifts on the 
cutoff parameter A that appears in the regulator function, Eq. ([2| . We vary A over a wide range, from 
0.5 GeV to 10 GeV (with n = 2 for the LO, NLO, and NNLO potentials, and n = 3 for N 3 LO). 

For partial-waves with short-range repulsion, convergence with increasing cutoff values is obtained 
without the use of any counterterm, i.e., convergence occurs "automatically". In fact, in these cases, 
any counterterm becomes ineffective for large cutoffs. Thus, no counterterm is used in partial-waves 
with short-range repulsion. 

For partial-waves with short-range attraction, one counterterm (contact term) is needed to ensure 
convergence. If we introduce a second counterterm per partial-wave, it turns out that this second 
parameter becomes ineffective for large cutoffs. Thus, we apply only one counterterm, which we use to 
fit the following empirical information. In S'-waves, we fit the scattering lengths; 1 Sq : a s — —23.740 fm, 
3 5i : at = 5.417 fm. In the other partial-waves with J < 2, we fit the phase shift at Ti ab = 50 MeV 
to the central value from the Nijmegen multi-energy np phase-shift analysis [5S]. For J > 3, we fit the 
phase-shift at 100 MeV or 200 MeV. 

In summary, at any order, either one or no counterterm is needed in each partial-wave for the 
infinite-cutoff renormalization. We show this in the left half of Table [T] for the various partial-wave 
states. The right half of Table [l] shows the number of counterterms according to Weinberg Counting. 
Obviously, there are large differences between the two schemes. The higher partial-wave states that 
are not listed in the table do not receive counterterms. 



9 



Table 1 Number of counterterms per partial-wave as required in two different renormalization schemes. 





Infinite-cutoff renormalization 


Weinberg Counting 


Partial- Wave 


LO 


NLO 


NNLO 


N^LO 


LO 


ATT /~\ /ATTVTT v*~\ 

NLO/NNLO 


N J LO 




f 


f 


f 


1 


1 


2 


4 


3 Po 


f 














1 


2 


L Pi 

















1 


2 


3 A 





1 


f 


1 





1 


2 


3 5x 


f 





f 


1 


1 


2 


4 


3 Si - 3 Di 








f 








1 


3 


3 £i 




















1 


L D 2 





f 


f 


1 








1 


3 D 2 


f 





f 











1 


3 P 2 


f 


1 


f 


1 





1 


2 


3 P 2 - 3 F 2 




















1 


3 F 2 























L F 3 























3 F 3 








f 


1 











3 D 3 








f 











1 


3 D 3 - 3 G 3 























3 G 3 























X G 4 








f 


1 











3 G 4 























3 F 4 








f 


1 












We have obtained convergence with increasing cutoff for all phase parameters with J < 6 in each 
order, NLO, NNLO, and N 3 LO. As evident from Table [lj in NNLO and N 3 LO, we need counterterms 
(contact terms) for the renormalization of the 3 F 3 , 3 i<4, and 1 G± waves. These terms are: 



(48) 
(49) 
(50) 



In Figures [T] and [2j we demonstrate the convergence of the phase parameters with increasing cutoff 
for the N 3 LO potential in partial- wave with J < 2. The black dotted curve is obtained for A = 0.5 GeV, 
the blue dashed curve for A = 1 GeV, the green dash-dot curve for A = 5 GeV, and the red solid curve 
for A = 10 GeV. The fact that the green and red curves are essentially indistinguishable in those 
figures, demonstrates that convergence for large cutoffs has occured. However, it is also clearly seen 
in these figures that, even though we are here at a relatively high order (N 3 LO), the cutoff-converged 
curves show large discrepancies with respect to the empirical phase shifts in several partial waves, 
particularly, 1 Sq, 3 Si, 3 -Di, 3 F>2, and £2- This issue will be further discussed in the next section. 



4 Order By Order Convergence 



Having accomplished the infinite-cutoff renormalization of the four orders we consider, it is now of 
interest to investigate the order-by-order convergence of these cutoff-converged cases. This is done in 
Figs- HE 




Fig. 1 Phase-shifts of neutron-proton scattering at order N 3 LO for total angular momentum J < 1 and 
laboratory kinetic energies below 250 MeV. The black dotted curve is obtained for A = 0.5 GeV, the blue 
dashed curve for A — 1 GeV, the green dash-dot curve for A = 5 GeV, and the red solid curve for A = 10 GeV. 
Note that the curves for A — 5 GeV and 10 GeV are, in general, indistinguishable on the scale of the figure. 
The filled and open circles represent the results fro m t he Nijmegan multi-energy np phase-shift analysis [58] 
and the VPI/GWU single-energy np analysis SM99 !59|, respectively. 



In Fig. [3j the S = 0, T = 1 np phase-shifts for L < 6 are displayed, where S denotes the total spin, 
T the total isospin, and L the orbital angular momentum of the two-nucleon system. Note that, in each 
order, the underlying analytic expression for the potential is the same except that it is decomposed 
into partial-waves with L = 0, 2, 4, and 6. The effect of this partial-wave decomposition is that the 
short-range part of the two-nucleon potential is increasingly suppressed with growing L; or in other 
words, the "centrifugal barrier" becomes larger with L. In the 1 S'o state, the phase-shifts are (almost) 
converged at N 3 LO, but do not reproduce the empirical phase-shifts. This is consistent with what was 
found in Ref. 42 ]. However, at the next higher S = 0, T = 1 partial- wave, the 1 D 2 , we observe that 
the phase-shift prediction converges to the empirical phase-shift values, and this is also true for all 
higher S = 0, T = 1 partial-waves shown in Fig. [3j 

In Fig. [4j we show the S = 0, T = state for various partial- waves up to L = 5. The situation is very 
similar to what we just discussed. In the lowest partial- wave, the 1 Pi, we have (almost) convergence, 
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Fig. 2 Same as Fig. [l] but J — 2 phase-shifts and J < 2 mixing parameters are shown 
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Fig. 3 Renormalized np phase-shifts at order LO (black dotted curve), NLO (blue dashed curve), NNLO 
(green dash-dotted curve), and N 3 LO (red solid curve). The S = 0, T = 1 phase shifts with L < 6 are shown 
for energies below 250 MeV. Filled and open circles are as described in Fig. [l] 
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Fig. 5 Renormalized S = 1, T = 1 uncoupled np phase-shifts with L < 5. Notation as in Fig. [3] 




though not to the empirical values. In higher partial- waves, however, the predictions converge to the 
experimental phase-shifts. This finishes the discussion of (spin) singlet states (5 = 0). 

Now we turn to (spin) triplet states (S = 1). Here, we need to distinguish between uncoupled 
and coupled partial-waves and will discuss the former first. The phase-shift predictions for uncoupled 
S = 1, T = 1 states are shown in Fig. [5] and for S = 1, T = in Fig. [6] Good convergence to the 
empirical phase-shifts is observed for all S — 1, T — 1 partial- waves (including the lowest one, 3 Pi). In 
contrast, for S = 1, T = 0, a strong divergence occurs at N 3 LO in the 3 £>2 state (first frame of Fig.[6|, 
while the higher S = 1, T = partial- waves converge well. 

Finally, we discuss the coupled cases. The T = 1 coupled partial-waves are displayed in Fig. [7} We 
include here the 3 Po state, such that we can compare it with its counterparts of larger L, namely 3 F 2 , 
3 i?4, and 3 Kq. It is clearly seen that the 3 Po shows no order-by-order convergence, while the 3 i<2 has 
near convergence and the 3 i?4 and 3 Kq states are fully converged to the empirical information. The 
associated coupled partial-waves show corresponding trends. 
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Fig. 7 The renormalized 3 Po and renormalized S — 1, T = 1 coupled np phase parameters for J < 6. Notation 
as in Fig. [3] 

The last set of partial- waves to be discussed are the S — 1, T — coupled states shown in Fig. [8j 
While there is no convergence with increasing orders for the J = 1 states, converged results and 
agreement with the experimental parameters are seen in the L = J + 1 partial-waves ( 3 G3 and 3 7s) 
and in the mixing parameters with J > 3 (i.e., £3 and £5). Contrary to this, the L = J — 1 waves 
( 3 Si, 3 £>3, 3 Gs) never converge. Even in the rather high partial-wave, 3 Gs, there is a large difference 
between the NNLO and N 3 LO predictions. This phenomenon may be related to the fact that the 3 G$ 
is the only higher partial- wave that disagrees with the empirical phase-shifts at N 3 LO, as was noticed 
already in Ref. |19j . At the present time, we do not have an explanation for this problem. 

The observations we have made in conjunction with Figs. [3][8] can be summarized as follows: Some 
lower partial-waves show no convergence and no order-by-order improvement. On the other hand, 
all higher partial-waves (except for the notorious 3 G5) do not only converge, they also converge to 
the empirical values. Note that in higher partial-waves the short-range part of the NN interaction 
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(equivalent to the high-momentum components of the interaction) is suppressed by the centrifugal 
barrier. This suggests that the long- and intermediate-range part of the NN potential up to N 3 LO is 
reasonable, while the short-range part may be, in part, un-physical. 

In the infinite-cutoff renormalization the potential is admitted up to unlimited momenta. However, 
the EFT this potential is derived from has validity only for momenta smaller than the chiral symmetry 
breaking scale A x » 1 GeV. The lack of order-by-order convergence and the discrepancies in lower 
partial-waves demonstrate that the potential should not be used beyond the limits of the effective 
theory (see Ref. [BU] for a related discussion). The conclusion then is that cutoffs should be limited to 



5 Summary and Conclusions 

We have investigated a particular scheme for the non-perturbative renormalization of the nucleon- 
nucleon (NN) potential based upon chiral effective field theory (chiral EFT). In the scheme applied, 
the cutoff parameter of the regulator is taken to infinity. 

Two vital requirements for constituting the legitimacy of an EFT are regulator independence and 
a power counting scheme that allows for order-by-order improvements of the predictions. We were 
able to achieve regulator independence (i.e., cutoff-independence) for all orders of chiral perturbartion 
theory considered (i.e., LO, NLO, NNLO, and N 3 LO) and all partial- waves up to J = 6. In general, 
cutoff-independence is seen for cutoff values above 5 GeV. 

However, in this investigation, we have also observed that large cutoffs impose limitations on the 
effectiveness of counterterms. In each partial- wave, either no counterterm (case of short-range repulsion) 
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or one counterterm (case of short-range attraction) is effective. Therefore, the power counting scheme 
implied by the infinite-cutoff method is considerably different from the one of Naive Dimensional 
Analysis or Weinberg Counting (cf. Table [IJ. As a consequence, order-by-order improvements of the 
predictions do not occur in several lower partial-waves. 

Thus, the chiral EFT approach to nuclear forces fails when renormalized by the infinite-cutoff 
method. This result may not come as a surprise considering that the EFT which we apply is designed 
for momenta below the chiral-symmetry breaking scale, A x w 1 GeV. Under the infinite-cutoff method, 
the potential contributes for momenta that are far beyond the hard scale of 1 GeV. Our results suggest 
that finite-cutoffs < f GeV should be used for the non-perturbative regularization of the chiral NN 
potential. With such cutoffs, all counterterms of Weinberg Counting are effective, and an order-by-order 
improvement of the predictions is to be expected |61J. 
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